Self-energy analysis of frequency-dependent conductivity: 
Application to Pb, Nb, and MgB 2 in normal state 
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We propose and demonstrate a microscopic way to analyze the frequency-dependent infrared 
conductivity: extraction of the electron self-energy from the inversion of experimentally measured 
infrared conductivity through the functional minimization and numerical iterations. The self-energy 
contains the full information on the coherent and incoherent parts of interacting electrons and, 
therefore, can describe their charge dynamics even when the quasi-particle concept is not valid. 
From the extracted self-energy, other physical properties such as the Raman intensity spectrum 
and the effective interaction between electrons can also be computed. We will first demonstrate 
that the self-energy analysis can be successfully implemented by fitting the frequency-dependent 
condcutivities of the simple metals such as Pb and Nb, and then calculating the effective interactions 
between electrons from the extracted self-energies and comparing them with those obtained from 
the tunneling experiments. We then present the self-energy analysis of the MgP)2 superconductors 
in normal state and clarify some of the controversies in their optical spectra. In particular, the 
small electron-phonon coupling constant obtained previously is attributed to an underestimate of 
the plasma frequency. 
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INTRODUCTION 

The frequency-dependent conductivity cr(u>) provides 
one of the most valuable and detailed information on the 
charge dynamics in a wide class of materials. It is ana- 
lyzed using either the one-component or two-component 
models EL The two-component model interprets the 
<j(lu) as arising from a combination of two types of car- 
riers, free and bound ones. The free carriers are mod- 
eled in terms of the Drude term and the bound ones in 
terms of various Lorentzian oscillators. However, the in- 
terpretation of the individual Lorentzian terms, most of 
which are usually due to inter-band contributions, is not 
straightforward. 

In the one-component picture, referred to as the ex- 
tended Drude model (EDM), on the other hand, the fre- 
quency dependence of the conductivity cr{u>) below inter- 
band contributions is described by extending the phe- 
nomenological parameters of the Drude model, the effec- 
tive mass m* and scattering rate 1/t, to be frequency 
dependent as Q 

ne 2 1 

cr(w) = 7TT7\ : 7r~Vi — ' W 

mi, 1/t{u>) — turn* (oj )/ nib 

where m& is the electron band mass. The EDM interprets 
the experimentally obtained complex conductivity <j{lu) 
in terms of 1/t(u>) and m*(u>) determined by 
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where lo p — (Anne 2 /m^) 1 / 2 is the plasma frequency, n 
the electron density, and e is the electron charge. The u p 
can be found by integrating the real part of the measured 
conductivity from the sum rule 



dw(j\{uj) = -lo 2 



(3) 
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where to max is the cutoff frequency above which inter- 
band contributions begin to contribute. The subscript 
1 and 2 refer to, respectively, the real and imaginary 
parts. The EDM has been successfully employed to ana- 
lyze <t(w) of the conventional metals as well as the heavy- 
fermions and high-T c cuprates [Q. 

For a class of correlated electron systems such as the 
ruthunates, Sr/CaRu03, however, the EDM breaks down 
and yields unphysical descriptions of the materials such 
as the negative effective mass || 0]. Similar behav- 
ior is also found for some molybdates, Sm2Mo20y and 
Nd2Mo2C>7 [||. These observations clearly signal the in- 
adequacy of the EDM and call for a new way of analyzing 
<r(u>) which can be applied to a wide class of materials. 
We, therefore, propose to analyze the frequency depen- 
dent infrared (IR) conductivity in terms of the electron 
self-energy £(u>) instead of the phenomenological param- 
eters of Eq. (||) . The electron self-energy contains the full 
information on the coherent and incoherent parts of inter- 
acting electrons and, therefore, can describe the charge 
dynamics even when the EDM or the Fermi liquid (FL) 
picture is no longer valid. In this self-energy analysis 
(SEA) method, the electron self-energy is extracted by 
inverting the experimentally measured infrared conduc- 
tivity through the functional minimization and numer- 
ical iterations as will be discussed in detail in Section 
III. The SEA may be considered as a microscopic gener- 
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alization of the EDM, which can be applied to analyze 
the frequency-dependent conductivity data of the non- 
Fermi liquids as well as the Fermi liquids. Even for the 
FL where the EDM is expected to work, the SEA can 
yield quantitatively more reliable results than the EDM, 
especially for the FL with a strong electron-phonon cou- 
pling. The present paper is mainly devoted to a detailed 
description of the SEA method and its applications to 
relatively simple metals. For more complicated cases of 
strongly correlated electron systems including cuprates, 
molybdates and ruthunates, by which we were originally 
motivated, we plan to report the SEA results separately 
elsewhere. For these systems, qualitatively different re- 
sults from the EDM results are expected. 

After the Introduction, we will discuss the frequency- 
dependent conductivity cr{u>) expressed in terms of the 
self-energy E(w) and its relation with the EDM in Sec. II. 
This will clarify the inherent limitations of the widely em- 
ployed EDM analysis for the frequency-dependent con- 
ductivity. We will then describe in Sec. Ill the formula- 
tion of the SEA method, which is reduced to the global 
minimization of a TV- variable function. The SEA method 
will be applied to experimental data of Pb, Nb and MgP>2 
and the results are presented in comparison with the 
EDM analysis in Sec. IV. In Sec. V, we give a brief sum- 
mary and some perspectives on the self-energy analysis 
method. 



OPTICAL CONDUCTIVITY AND SELF-ENERGY 



where E 2 mp is the frequency-independent contribution 
from impurities, /^(w) is the imaginary part of the effec- 
tive interaction satisfying P 2 {—<jj) = —P 2 (uj) and T is the 
temperature. The isotropically weighted phonon density 
of states for electron-phonon coupled systems is given by 



a 2 F(uj) = -P 2 (u). 

7T 



(6) 



Once the self-energy is known, we use Eq. ([s]) to find the 
effective interaction spectrum J^C^O using a derivative 
with respect to lo at low T or a convolution after Fourier 
transformations [Q. 

In the zero temperature limit T = 0, the ct{uj) of Eq. 
(0) is reduced to 
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to — — e) — E(e) ' 



(7) 



which means that in the low T limit, only the self-energy 
S(e) between < e < lu contributes to the conductivity 
at to. This may be interpreted as <t{uj) being an "average" 
of — E(e) — T,(lo — e)] between and uj. This, in turn, 
suggests that the information from the EDM analysis, 
which is directly obtained from a(ui) using Eq. (||), is 
an average of the corresponding quantity from the SEA. 
This will be discussed in more detail below. The EDM 
can be obtained from Eq. (||) in an appropriate limit. 
The a(oj) of Eq. (|) is reduced to the EDM form of Eq. 
(Q) with m*(uj)/mb = 1 + A(u) provided that 



Yi{u> — e) + S(e) « — lo\(lu) — i/riui) 



(8) 



The frequency-dependent conductivity <j(uj) can be ob- 
tained from the current-current correlation function, and 
is written in terms of the electron self-energy E(oj) as 



a[u) = — - — 

47TW 



f(e-u)-f(e) 



^ d6 cj - S(w - e) 
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where /(e) is the Fermi distribution function. We as- 
sumed a constant density of states over an infinite band- 
width and no long-range order. It is also assumed that 
the momentum dependence is much weaker compared 
with the frequency dependence, S(fc,w) = as in 

the dynamical mean-field theory, which renders the ver- 
tex correction vanish in the current-current correlation 
function ||. 

The electron self-energy represents the effects of elec- 
tron interaction with various excitations in a system. The 
imaginary part of the self-energy can be written as 



EaM = - / 

J — ( 



coth[- 



tanh 



2T 



p 2 m 
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For a FL, where Ei(w) w — \(lu)u> and 



E™' — Tuj 2 , this condition can be satisfied for 



is satisfied 
E 2 (w) 

small e and u>, if A(w) has a weak w-dependence. There- 
fore, the EDM can give a satisfactory description of a(ui] 
for weak-coupling FL, where 



1 ~-2£ 2 (u,), m * {u) 



r( W ) 



mb 



1 + A(w) 



In the weak-coupling limit, where Eq. (g) is well satisfied, 
the optical scattering rate 1/t(lo) can be approximately 
written as 10, 0. ml 
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where 1/Ti mp is the impurity contribution and a 2 r F(uj) is 
a phonon density of states weighted by the amplitude for 
large-angle scattering on the Fermi surface, which has the 
same spectral structure as a 2 F{u>), but their amplitudes 
can be lower. In this paper we will not distinguish the 
a 2 F(u>) and a 2 r F(oj) from now on. We note that there 
were several previous attempts to invert a 2 F{uj) H, fflj. 
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For instance, one may obtain a 2 F(uj) in the T = limit 
of Eq. (jlC|) using 



2nduj 2 



1 ^ d 2 
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cr(w) 
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where the second expression follows by using the EDM 
of Eq. (|l|). This formula, however, has limitations to be 
applied to the experimental IR data because large error 
bars are inevitable from the double differentials and it is 
valid only when EDM is valid and at T = @ . 

For strong-coupling FL where A(w) has a significant 
^-dependence, analysis based on the EDM becomes less 
accurate. For the marginal Fermi liquid, where £i(cj) ~ 
win and £2(0;) ~ ~M, the EDM is expected to give 
somewhat less reliable description because £i(w), unlike 
FL, deviates from the linearity in u>. The situation be- 
comes worse for non- Fermi liquid, where £i(w) ~ — ui 1 ~ a 



and SMw) 



(0 < a < 1), and the EDM may 



give misleading and qualitatively incorrect descriptions 
as were observed in the ruthnates H H • 



FORMULATION OF SELF-ENERGY ANALYSIS 

We now present how one can analyze the frequency- 
dependent infrared conductivity with the formula of Eq. 
(||). As is explained below, the problem is reduced to a 
global minimization of a N- variable function. Using Eq. 
(Q), it is simple to calculate the conductivity cr(uj) from 
a given self-energy £(w). What we are trying to do here 
is exactly the inverse of that: We wish to extract £(w) 
from an experimentally measured <r(w). 

The SEA is implemented by defining the functional 
W£ 2 ] as 



W\T, 2 



dui [ui (w) 



exp / \u 



(12) 



where a^ xp and o~\ are the real parts of the experimental 
data and calculated conductivity using Eq. (Q), respec- 
tively. Since the real part of self-energy £i(w) can be 
obtained from the imaginary part £2(0;) by the Kramers- 
Kronig transformation, we can consider £ 2 (w) as the 
only independent function. The functional W is posi- 
tive definite, and has the global minimum of zero for the 
self-energy which reproduces the experimental conduc- 
tivity data. The self-energies defined at the N discrete 
frequencies, Xi = £2(^1) f° r i — 1, * • • , iV (wn = w c ), 
are taken as the N independent variables of the func- 
tion W(x\,X2, • ■ • ,xn) = \V\S2]. The cutoff frequency 
u! c is not necessarily equal to the io max of Eq. (|J). Now, 
the problem is reduced to a global minimization of a N- 
variable function. 

We note that, depending on the problems, other forms 
of the functional W may yield better results. For good 



metals, for instance, whose mid-IR conductivity is very 
small compared with far-IR region, 



W = 



duj [ln<7i(u;) — In a1 xp (lu)] 



(13) 



works better than the form of Eq. (JT^) . For the region 
where to > u c , which we need to know for the Kramers- 
Kronig transformation, £2(w) is taken as constant. A 
constant £2^) for to > oj c corresponds to the reflectiv- 
ity R(uS) ~ 1/w 4 , which is consistent with the standard 
procedure in the IR experiments. 

The global minimization of W is achieved via the func- 
tional derivative and numerical iterations. We start with 
an initial configuration of x^\ A good initial guess can 



be 



,(°) - 



obtained from EDM analysis, or a neg- 



ative constant value for the whole frequency range. Then, 
we move to new Xi along the steepest descent direction 
of W. 



(new) (old) 
X) — X) — S 



dW 

dxj 



(14) 



where represents the functional derivative of the 
functional ^[£2] with respect to T,2(u>) given by 

dW 



dxi 



<5eri(w') 1 



de *^y) p _L (15) 



^£2(0;) 7r J_ 00 5£i(e) uj — e 

where V stands for the principal value, and the use is 
made of f3M £ ' 1 , = -V—^—. The step size s is chosen such 
that the W is maximally decreased along the steepest de- 
scent direction. The cj 2 , which sets the scale of a(ui) of 
Eq. (Qj is updated at each iteration such that the calcu- 
lated spectral sum is equal to the experimental spectral 
sum up to uj c . Note that this way of determining oj p does 
not require that uj c — Lo max of Eq. (|^) and, therefore, the 
extracted uj p is almost independent of io max and more 
reliable than that from the sum rule of Eq. (§). This en- 
ables us to determine lu p more systematically as we will 
discuss in more detail below. 

In general, the global minimization of a few hundred 
independent variables is an extremely demanding prob- 
lem. In the present case, however, it is rendered tractable 
because of the following observations: (1) We know the 
value of the global minimum unlike general global mini- 
mization problems. It is exactly zero. (2) We have some 
ideas about the physically meaningful form of the £2(w). 
It should be a continuous function of the frequency and 
negative definite. (3) We have a better way to escape 
from local minima than trying a new random starting 
point, xf: Take£ 2 neuj) M-£2 OW) H oc oi(w)-of (w). 
An estimation of this process can be obtained in follow- 
ing way. The contribution to the conductivity at a given 
frequency lj is dominated by the self-energy below u> in 
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the low temperature limit due to the thermal factor of 
[f(e -u)- f(e)]/u of Eq. @ as discussed in Eq. (@). 
By modifying the self-energy to be proportional to the 
difference between calculated conductivity and the ex- 
perimental data at each frequency, we can continue our 
minimization modifying the 'wrong' region without spoil- 
ing a 'good' region. 

From the extracted other physical properties 

such as the plasma frequency, effective interaction be- 
tween electrons, Raman spectra, and inelastic neutron 
scattering intensity can also be calculated. For instance, 
we use Eq. (||) to find the effective interaction spectrum 
P2(u) from an extracted ^(u) which was discussed in 
the previous section. An important byproduct of the 
SEA is that the plasma frequency u p can be determined 
more accurately than has been hitherto practiced, which 
is necessary to determine the Drude parameters of Eq. 
(^) . Conventionally, u p is determined by uj max of Eq. (|^) 
which, in turn, is taken such that the sum as a func- 
tion of u max has the smallest slope. This procedure 
can be problematic especially for the materials without 
a sharp plasma edge in the reflectivity data. In con- 
trast, SEA works in this case as well because the deter- 
mination of u p , which sets the scale of cr(u), is almost 
independent of the u c of the functional W. For instance, 
one can determine u p with only very restricted data be- 
tween < u < u c (<C u max ). This will be illustrated for 
Nb below. 

Another byproduct closely related to the u p determi- 
nation is that we may separate the intra- and inter-band 
contributions more systematically. We can extract the 
self-energy of effective single-component carriers by fit- 
ting the low-frequency experimental conductivity up to 
u c w u max . The conductivity a(u) calculated by sub- 
stituting the extracted £(u>) into Eq. (||) is intra-band 
conductivity. This process, however, is more sensitive to 
u c than the u p determination. 

We will demonstrate in the next section that the SEA 
is straightforward to implement to extract the self-energy 
from experimental IR conductivity without the local min- 
imum problem. We have found that the SEA yields 
the same solution £ 2 (w) from almost any initial config- 
uration. The obtained solution T,2(u), therefore, seems 
unique. 

APPLICATIONS TO EXPERIMENTAL DATA 

We will now apply the SEA developed in the previous 
section to real materials, and demonstrate that it can be 
readily employed to analyze the frequency-dependent IR 
experimental data. In the first two parts of this section, 
we will analyze the simple metals, Pb and Nb. From 
the SEA, we can extract the self-energy of the materi- 
als, and from the extracted self-energy we can find the 
ct 2 F(u) using Eqs. (||) and (||). It can be compared with 



the measured a 2 F(u) well established from the tunneling 
experiments jl^] . But for such a good metal, it is a very 
demanding task to measure the conductivity in far-IR re- 
gion, since the reflectivity is very close to 100 % there. 
Thus there exist few published conductivity data for Pb 
and Nb. For Pb, we take the experimental a 2 F(u) ob- 
tained from the tunneling experiment and calculate the 
self-energy from Eqs. (0) and (||), and then calculate the 
conductivity from Eq. (0). This is taken as the "experi- 
mental" conductivity a exp (u) for Pb. After that, the real 
experimental data of Nb is analyzed jl3| . The results are 
compared with those obtained from the tunneling spec- 
tra. The last part of this section is devoted to the SEA 
for the normal state IR conductivity data of MgB 2 super- 
conductor. We will argue from the SEA results that the 
small electron-phonon coupling constant A extracted pre- 
viously from the T-dependence of the resistivity p(T) and 
EDM analysis of cr(u), which is too small to account for 
the superconducting transition temperature T c , is most 
likely due to an underestimate of the plasma frequency 
u p . 

Application to Pb: Generated data 

We take Pb as a test case of the proposed SEA be- 
cause its a 2 F(u) is well established from the tunneling 
p2| . The electron-phonon coupling constant A of this 
metal is known to be w 1.5. The frequency-dependent 
IR conductivity of Pb, however, is not available. There- 
fore, the IR conductivity data was generated with Eqs. 
(U) and (||) using the tunneling a 2 F(u) at 200 frequen- 
cies up to u c = 15 meV with an equal spacing. We took 
the self-energy due to impurity scattering £ 2 mp = — 1 
meV and T = 3 K as representative values. The gener- 
ated conductivity in this way, shown in the left column 
of Fig. 1 with open circle, is taken as the "experimen- 
tal" data al xp (u). It was then fitted by the SEA as ex- 
plained above to extract the self-energy T, 2 (ui) without, 
of course, any information about the self-energy used to 
generate the conductivity. The results are shown in the 
solid line in the left column of Fig. 1. The extracted and 
experimental conductivities are almost indistinguishable. 
Note that the width (HWHM) of ai(u) is not given by 
— 2^2(01 = 0) because of the substantial frequency de- 
pendence of £2(0;). 

In the right column of Fig. 1, we show the real and 
imaginary parts of the extracted self-energy in the forms 
of — 2S 2 (^) and 1 — Y,i(u)/u, along with the generated 
experimental data. The extracted and experimental data 
are shown, respectively, by the solid line and open cir- 
cles. They are practically indistinguishable. For com- 
parison, we also plotted the EDM analysis of the ex- 
perimental conductivity data, 1/t{u) and m*{u)/m , in 
dotted lines, using Eq. (||). As discussed in Eq. (^), 
1/t{u) andm*(w)/m(, correspond, respectively, —2Y,2(u) 
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and 1 — Ei(ti))/aj. The right panel demonstrate, how- 
ever, that they can substantially deviate from each other 
even for the relatively simple metals. The EDM analy- 
sis in general yield smoother frequency dependences for 
1/t(u>) and m*(u))/mt > . This is expected because 1/t(u>) 
is an "average" of the £2(0;) as discussed previously in 
Section II. The curve shown by the thin solid line in the 
lower picture was calculated by using the approximate 



formula Eq. (10), which is in a good agreement with the 
1/t(w) from the EDM analysis @. 

From the extracted self-energy, a 2 F(cu) can be ob- 
tained from a derivative with respect to u>, a 2 F(uj) — 
— \ 9Ti q^ , which is valid at low T (dotted curve) or a 
convolution of Eq. (||) (thin solid curve) JtJ as shown in 
the upper left inset of Fig. 1. The experimental a 2 F(ui) is 
shown in the thick solid line. The extracted a 2 F(u>) from 
the convolution is again almost indistinguishable from 
the experimental one. These results from the SEA show 
substantial improvements over the previous attempts Q 
based on the assumptions of both T = and weak- 
coupling limit. We argue from this example that even 
for the simple metals, for which the EDM can give the 
qualitatively valid description, the SEA can provide more 
accurate and reliable description of the material, which 
can be quite different from the EDM results. 



Application to Nb 

We now proceed to apply the present method to the 
experimental IR data of Nb. We used the far-IR conduc- 
tivity data measured by Pronin et al. in the normal state 
at 9 K Their far-IR data are available from 84 up 
to 300 cm -1 (with constant data in the DC limit mea- 
sured via a different method), which is a good example 
to demonstrate that the oj p can be obtained without a 
full intra-band spectrum in the SEA. 

We took lo c = 250 cm" 1 and N = 200. For < to < 
84 meV, we constrain ^(w) to remain constant. This 
implies vanishing a 2 F(u>) in the region as shown in 
the upper right inset of Fig. 2 and, consequently, a re- 
duced A, and the Drude behavior in eri(w). The fitted 
(solid curve) and experimental (crosses) <7i(w) are plot- 
ted in Fig. 2. The electron-phonon coupling constant 
A = —dT,i(u;)/duj\ U j^o from the SEA is 0.51 which is 
somewhat reduced compared with the experimental value 
of 0.9 — 1 Q as discussed above. The calculated plasma 
frequency is 7.8 eV while the experimental value is 7.2 
cV [|l3|. In the lower left inset, the extracted £2(0;) is 
shown. The above results clearly demonstrate that the 
SEA works also for the case with only restricted data. 

These two examples of the SEA establish that the 
method can indeed be applied to analyze the frequency 
conductivity and it can provide the most microscopic in- 
formation of interacting electron systems, the self-energy 
£(w) and the effective interaction P(oj). The SEA can 



provide for the Fermi liquids more reliable and accurate 
information than the conventional EDM analysis. For 
the non-Fermi liquids, it is expected to provide qualita- 
tively different information which is not accessible with 
the EDM. 



Application to MgB2 

Let us now analyze the IR data of normal state c-axis 
oriented MgP>2 film measured by Tu et al. E3, which 
we regard as being due to intra-band excitations of an 
effective single band system 0. In Fig. 3(a), the experi- 
mental conductivities of Tu et al. at T = 45 K and 295 
K are shown together with the fitted conductivities us- 
ing the SEA with N = 300 and to c = 6000 cm" 1 . The 
solid curves represent the fitted conductivities, and the 
dashed and dotted curves, respectively, the experimental 
ones at T = 295 K and T — 45 K. The small discrepan- 
cies are due to phonons. In Fig. 3(b), the results from 
the EDM analysis are shown in the left column, and those 
from the present SEA in the right column. The solid and 
dashed curves are, respectively, for T — 45 K and 295 
K. Tu et al. found from the sum rule of Eq. (||) that u p 
is 14750 cm -1 which yields, through the EDM, a very 
weak electron-phonon coupling as shown in the left col- 
umn of Fig. 3(b). They obtained Xt r ~ 0.13 using the 
Bloch-Griineisen formula to analyze the T-dependence of 
the resistivity, which is consistent with the EDM analysis 
p4] , To the extent that \ tr ~ A, it seems too small to 
account for the superconducting transition temperature 
T c = 39 K. On the other hand, we found from the SEA 
that A w 0.56 and 0.41 at, respectively, T — 45 and 295 
K as shown in Fig. 3(b), which are substantially larger 
than what Tu et al. found. The local density approxima- 
tion calculation yields Aj r « 0.6 [16). One way of seeing 
the discrepancy between EDM and SEA is that the uj p 
from SEA is larger than that from EDM. We found that 
lu p = 16690 cm" 1 at T = 45 K, and 16740 cm" 1 at 
T = 295 K. If the enhanced u> p are used in the EDM, the 



resulting X tr are in good agreement with the SEA. As far 
as the EDM and the T-dependence of the resistivity are 
concerned, the enhanced uj p resolves the problem of the 
small X tr . 

With the SEA, we may go further and perform the 
spectral analysis to see the frequency range that con- 
tributes to A, which can not be carried out with the EDM 
analysis. A = 2 f™ dtl "*^ , where a 2 F{uj) = ±P 2 {uj) 
may be obtained from the extracted £2^) using Eq. (|^). 
The — 2£2(w) is shown in the right column of Fig. 3(b). 
The extracted a 2 F{u>) is shown in Fig. 4 together with 
those from tunneling Jl7j (dashed line) and LDA calcula- 
tion [ fl6| (thin solid line). They are shown in a wider 
frequency range in the inset. The extracted a 2 F(w) 
from SEA are characterized by two frequency regions 
which make dominant contributions to A; around w w 
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70 and 300 meV, which is, interestingly, consistent with 
the model proposed by Marsiglio Q . The low frequency 
region around 70 meV is the contribution from the E 2g 
phonon mode, while the nature of the high frequency 
region is not clear. The two regions contribute almost 
equally to A: Xphonon ~ 0.36 (0.21) and the total A is 
0.56 (0.41) for T = 45 K (T = 295 K). 

The extracted a 2 F (w), compared with the local den- 
sity approximation (LDA) calculation, correctly captured 
the main E 2g contribution but missed smaller contribu- 
tions from other phonon modes; among the modes in 
the a 2 F(uj) obtained from the LDA calculation shown in 
Fig. 4, the modes whose a 2 F{uj) are smaller than ~ 0.5 
are absent in the SEA results. This, we suspect, may 
be a consequence of the broad features around 160 and 
880 cm -1 in the experimental <r(w), which were not pre- 
dicted by phonon calculations Q and possibly due to 
MgO impurities. We did not eliminate these contribu- 
tions for the present calculations. This may smear out 
otherwise sharper frequency dependence of o~{oj), and re- 
duce the electron- phonon coupling constant. 

Apart from these discrepances, the SEA successively 
describes the frequency-dependent conductivity of the 
MgB 2 and yields an increased A than the previous es- 
timate. The SEA suggests, which made more accurate 
determination of w p possible, that the small A reported 
in the c-axis oriented sample is most likely due to under- 
estimated uip, and the total A from the SEA is substan- 
tially larger. However, the Xphonon ~ A/2 still seems a 
bit too small to account for the T c . 

SUMMARY AND DISCUSSIONS 

In this paper we have proposed and demonstrated a 
microscopic way to analyze the frequency-dependent in- 
frared conductivity, referred to as self-energy analysis, 
which is valid for the non-Fermim liquids as well as the 
Fermi liquids. Additional advantage of the self-energy 
analysis is that the plasma frequency uj p can be obtained 
with a better accuracy. Through the self-energy analysis, 
we extracted the electron self-energy from the inversion of 
experimentally measured infrared conductivity, and the 
effective interaction between electrons from the extracted 
self-energy. After we demonstrated that the self-energy 
analysis method can be successfully applied for simple 
metals like Pb and Nb, we applied the method to the IR 
conductivity data of normal state MgB2. We have found 
that the A from the self-energy analysis is substantially 
larger than that obtained from the conventional analy- 
sis of the T-dependence of the resistivity and extended 
Drude model. The discrepancies between the self-energy 
analysis and conventional method for MgB2 is attributed 
to the underestimated u> p from the conventional method. 
However, the Xphonon from the self-energy analysis still 
seems a bit small to account for the T c . 



Now that we have demonstrated that the self-energy 
analysis of the frequency-dependent infrared conductiv- 
ity really works and can be very powerful, some con- 
cluding remarks and outlooks are in order. First, one 
may wonder if the solution to the global minimization 
of the functional H^[S 2 ] is unique. The answer seems 
positive: the converged solutions, with different initial 
configurations, of given conductivity data all agree with 
each other. This means that the extracted self-energy 
S2(w) is indeed unique. Second, we have also checked if 
there is any spurious feature from the Kramers-Kronig 
transformation used in the present work to obtain the 
real part of the self-energy from the imaginary part. We 
therefore have fitted the complex conductivity by treat- 
ing both the real and imaginary parts of S(w) as inde- 
pendent variables, which eliminates the use of Kramers- 
Kronig. This gives the same £2(^0 with the proce- 
dure using iy[E 2 ], demonstrating the reliability of the 
present method. Third, we plan to report results of 
the self-energy analysis of the conductivity for other cor- 
related electron systems such as the high T c supercon- 
ductors, Sni2/Nd2Mo207, and Ca/SrRuC>3 mentioned in 
the introduction. The differences between the self-energy 
analysis and the extended Drude model analysis are ex- 
pected to be quantitative for the Fermi liquids, but be 
qualitative for the non-Fermi liquids. It will therefore be 
very interesting to see what information the self-energy 
analysis provide for the strongly correlated electron sys- 
tems. Also among the plan are extending the self-energy 
analysis to non constant density of states and broken 
symmetry states, and doing the analysis for other two- 
particle probes such as the electronic Raman and inelas- 
tic neutron spectra. We believe that the self-energy anal- 
ysis seems timely and urgent in view of the mounting 
interests in the correlated electron or non-Fermi liquid 
systems for which the phenomenological analyses may 
yield inadequate results. 
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FIG. 1: Left: cn(w) is plotted for Pb with E l 2 mp = -1 meV 
and at T — 3 K. The extracted and generated conductivities 
are shown, respectively, in solid line and open circles. The 
dashed and thin solid curves in the inset were obtained from, 
respectively, a 2 F(ui) = ~~ ~~§£r^ an d the convolution using 
Eq. (^), and the latter is almost indistinguishable with the 
input tunneling data (solid line). Right: The 1 — Ei(w)/w and 
— 2£2(cj) obtained from the SEA are plotted, respectively, in 
the upper and lower panels. The extracted (generated) ones 
are shown in solid line (open circles). These SEA analysis 
results are compared with the EDM results. The dotted lines 
represents 1/t(u) and m*(uj)/mt,, respectively, obtained from 
EDM using the generated conductivity. 1/t(u) (thin solid 
line) calculated by using Eq. (|Io|) is also shown in the lower 
right panel for comparison. 
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FIG. 2: The calculated and experimental cri(a>) of Nb are 
shown with solid curve and crosses, respectively. In the insets, 
the extracted o?F{ui) and E2(oj) are shown along with the 
tunneling o?F{ui) (dashed). 




FIG. 3: (a) The calculated and measured ctl(w) of MgB2 
for T = 45 and 295 K. The small discrepancies are due to 
phonons. (b) Left: extended Drude analysis of experimental 
data indicating the weak electron- phonon coupling. Right: 
the corresponding quantities from the self-energy analysis. 
Note that m* /nib and l/r(u) show much smoother frequency 
dependence compared with the corresponding 1 — Si (a;) and 
— 2E2(w) as noted previously for Pb in Fig. 1. 




FIG. 4: The extracted a 2 F(cu) from E 2 (w) of Fig. 3 along with 
those from LDA and tunneling. The extracted a 2 F(u>) shown 
in a wider frequency range in the inset are characterized by 
two dominant contributions around u> £s 70 and 300 meV. 



